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Abstract 

We use a tunable laser ARPES to study the electronic properties of the prototypical multiband 
BCS superconductor MgEG Our data reveal a strong renormalization of the dispersion (kink) 
at ~65 meV, which is caused by coupling of electrons to the Y>g phonon mode. In contrast to 
cuprates, the 65 meV kink in MgB 2 does not change significantly across T c . More interestingly, 
we observe strong coupling to a second, lower energy collective mode at binding energy of 10 meV. 
This excitation vanishes above T c and is likely a signature of the elusive Leggett mode. 

PACS numbers: 74.25.Jb, 74.72.Hs, 79.60.Bm 
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In conventional superconductors, the pairing is mediated by phonons and favored by 
strong electron-phonon coupling, as described by Bardeen-Cooper-Schrieffer (BCS) theorjf 1 . 
This strong electron-phonon coupling in general gives rise to a renormalization of the band 
dispersion called a “kink” and an abrupt change of quasiparticle lifetime at an energy re¬ 
lated to the phonon frequency, fh This idea has been extended to unconventional super¬ 
conductors, where the mechanism of pairing is unknown, and the coupling of electrons to 
several collective excitations was reported 4 ' 9 . Their origin and relation to pairing is still 
debated. For example the ”70meV” kink in Bi2212 cuprate is strongest at the antinode 
and vanishes above T.P a behavior that resembles the magnetic resonance mode reported 
by inelastic neutron scattering 10 ; however, the electron-phonon interaction may have similar 
characteristics 1112 . On the other hand, the kinks in the dispersion along the nodal direction 
in single layer Bi2201 do not seem to change significantly with temperature 14 . Surprisingly, 
there is little data on dispersion renormalization effects in conventional superconductors; 
several low energy kinks have been reported in NbSe^, but the situation is complicated 
by the presence of a charge density wave phase coexisting with superconductivity, and the 
measurements were carried out only at low temperatures. In other materials, difficulties in 
clearly observing three dimensional band dispersions and low transition temperatures are 
limiting factors. MgB 2 15 20 is a notable exception: it is a layered material with multi-gap, 
phonon-mediated superconductivity at T c = 39K with some quasi-two-dimensional bands, 
making it an ideal candidate to study the temperature dependence of the dispersion renor¬ 
malization due to electron-phonon coupling. LDA calculations 17 predict four bands crossing 
Fermi level in MgB 2 : two quasi 2D cr-bands from p x , p y orbitals around T and two 3D 7 r- 
bands from p z orbital 2731 32 . The superconductivity is believed to caused by the E 2g phonon 
mode at 75meV that couples strongly to the 2D cr-bands, but more weakly to the vr-bands 1 ^, 
leading to two different gaps, A a = 6.5 meV and A n = 1.5meV, as revealed by previous 
tunneling 26 and ARPES studied. Inelastic neutron scattering studies have reported optical 
phonon modes at ~35, 55, 75, 85 and 100 meV 28 , while Raman scattering reports a single, 
broad asymmetric peak at 75meV, attributed to the E 2g mode, as well as a sharp peak 
at 2A a = 12meV due to a pair breaking excitatiorP 9 . As a bonus, due to its multi-gap 
nature, MgB 2 also contains another exotic collective mode that can couple to the electrons: 
the Leggett modP 1 25 . This mode is a longitudinal fluctuation corresponding to equal and 
opposite displacements of the two condensates. As it is a neutral excitation, it is not pushed 
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up to the plasma frequency like the Bogoliubov-Anderson mode, however, it is predicted to 
have a mass in-between 2A 7r and 2A a for MgB 2 , where it will be partially damped by de¬ 
caying into 7T quasiparticlcs- 0 '. As the superconductivity is believed to be phonon mediated 
and due to intraband pairing, the two gaps are expected to have the same relative sign, 
giving rise to s ++ pairing, by contrast to the iron based superconductors, where interband 
repulsive pairing likely leads to s± pairing and a Leggett mode is not expected. 

High resolution band dispersion data demonstrating the coupling of the conduction elec¬ 
trons to collective excitations such as phonons or Leggett modes are not available in the 
published literature. One of the main reasons is that high quality single crystals need to 
be synthesized under high pressure, resulting in rather small <500 pm crystals, which are 
difficult to measure in traditional ARPES setups. We use tunable laser ARPES to study the 
electronic properties of MgB 2 multiband superconductivity. The use of low photon energy 
increases significantly the bulk sensitivity due to increased escape depth and momentum 
resolution due to increased A -1 /deg ratio. The ability to focus the laser beam down to 
~30pm enables measurement of very small single crystals and also helps to improve the 
momentum resolution. We find evidence for strong coupling of conduction electrons to a 75 
meV acoustic phonon with A estimated at ~ 1.3 that persists above T c , unlike those in the 
cuprates. In fact we observe no significant changes with temperature up to 65K, more than 
50% above T c . Furthermore, we observe a very strong renormalization of the dispersion 
in the superconducting state at ~ 10 meV. Instead of the expected Bogoliubov-like back 
bending of the dispersion, a sharp, non-dispersive quasiparticle peak centered at the gap 
energy of 6.5 meV is present and is separated by a dip from the high energy spectral weight. 
All these features vanish above T c as expected for a superfluid excitation and are likely due 
to interaction of electrons with a Leggett mode. 

MgB 2 single crystals with T c = 39 K and typical size of 0.5 x 0.5 x 0.3 mm 3 were grown 
in Ames Laboratory by a high pressure synthesis technique similar to that outlined in Ref. 33 
using pure n B isotope. Optimally doped Bi 2 Sr 2 CaCu 2 0s+5 (Bi2212) single crystals with 
T C =93K (OP93K) were grown by the conventional floating-zone (FZ) technique. Sample 
were cleaved in situ at a base pressure lower than 8 x 10 -11 Torr. ARPES measurements 
were carried out using a laboratory-based system consisting of a Scienta R8000 electron 
analyzer and tunable VUV laser light soured 4 . All data were acquired using a photon energy 
of 6.7 eV. The energy resolution of the analyzer was set at 1 meV and angular resolution was 
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0.13° and ~ 0.5° along and perpendicular to the direction of the analyzer slits, respectively. 
Samples were cooled using a closed cycle He refrigerator and temperature was measured 
using a silicon-diode sensor mounted on the sample holder. The energy corresponding to 
the chemical potential was determined from the Fermi edge of a polycrystalline Au reference 
in electrical contact with the sample. The aging effect was checked by recycle measurements. 
The consistency of the data was confirmed by measuring several samples. 

The Fermi surface data are shown in Fig. la, with a schematic plot of the Brillouin zone 
inset. We used the peak position of MDCs at E F to quantitatively extract k^. Results are 
superposed as red points on image data. Both a FS sheets are round with |kp| ~ 0.2 A -1 
and ~ 0.25 A -1 respectively. If we ignore the small warping of these two a sheets along 
k 2 , the contribution to carrier concentration would be 0.069 holes for inner FS sheet and 
0.108 holes for outer one. The measured area of the outer FS is consistent with previous 
quantum oscillation results, while the inner one is slightly larger 36 . Fig. 1 (b, c) show 
ARPES intensity plots at two different cuts in the Brillouin zone. In the data measured 
along cut # 1 (Fig. lb) the Fermi crossing for both a sheets are clearly visible, whereas in the 
cut #2, along symmetry axis (Fig. lc), intensity of the inner cr-band is strongly suppressed 
due to the matrix element effect. In order to quantitatively analyze the renormalization 
effects due to the collective modes, we fit the MDCs of each data set with Lorenzians and 
plot such extracted dispersion as blue and red lines in Fig. 1 (b) and (c). It should be 
noted that MDC peaks do not reflect the dispersion at very low energies in the presence of 
the superconducting gap therefore the fitting is carried out only for binding energies larger 
than ~2A. We estimate bare dispersion by extrapolation from higher binding energies and 
plot them as dashed lines. In all three data sets, a very pronounced kink structure is clearly 
visible (indicated by arrows), where the renormalized dispersion deviates from bare estimate. 
In this case, the renormalization of the dispersion (kink) is the fingerprint of coupling the 
conduction electrons to phonon mode(s). This coupling is rather strong and the energy 
distribution curves (EDCs) at k F develop a dip (right side of panel 1 (c)). The kink in 
dispersion and dip in EDC occurs at ~ -(70-75) meV. By comparison with Raman datsP^, 
we can conclude that it is due to the E 2g optical phonon responsible for pairing in MgB 2 19 . 
We subtract an estimated bare dispersion from the measured one to obtain the approximate 
E'(u;), and use the widths of the Lorentzian MDC fits to obtain E "(u). These quantities 
are plotted in panel 1 (d). Not surprisingly, E'(cj) has a peak at -70 meV, which is very 
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FIG. 1. (Color online) Renormalization effects (kink) in cr-band. (a) Measured FS at 40K. Blue and 
red circles mark the Fermi momentum extracted from MDC peak position at E^. Two red arrows 
mark locations of measured cuts in the momentum space. Insert is a schematic diagram of the FS 
and Brillouin zone, (b) ARPES intensity along cut #1 at 15K. Solid lines are dispersions obtained 
by MDC fits for the two bands. Black dashed lines signify bare band dispersion extrapolated from 
higher binding energies, (c) left: same as (b) but along cut #2; right: EDC at for outer cr-band. 
(d) Energy dependence of the MDC width (blue solid line) and effective real part of self-energy 
(red dotted line) for data in panel (c). Arrows mark the energy location of the kink and associated 
features. 


5 






















close to the phonon energy. E"(c<j) rises rapidly with binding energy and has a mid-point of 
the step roughly located at the same energy as the peak in £'(u;). The line shape of both 
curves is consistent with them being related by Kramers-Kronig. This is further verified by 
simulations presented in Supplemental Material. 

The dispersion data extracted by fitting the positions of the MDC peaks for several 
temperatures both below and above T c are shown in Fig. 2(a). There are no significant 
changes to the kink structure across T c as evident from overlay of the high temperature curve 
(dashed red line) onto the lowest temperature one. The only noticeable changes occur at low 
energies due to the opening of the superconducting gap. Following the procedure outlined 
above, we extract £'(u;) for each temperature, and plot these in Fig. 2(b). The peak 
position does move to lower energies as the temperature decreases, starting at T c , consistent 
with the expectation that the peak frequency increases from 0 —» Q + A. However, the 
magnitude of the shift, 3meV is significantly smaller than the expected A^ = 6.5meV. In 
fact, this shift should be even larger as the screening electrons are gapped out and E 2g 
phonon hardens below T c . This hardening was predicted to be ~ 10% or 7meV^, but 
Raman measurements find that it only shifts by AQ « 2.5meV 37 . So the overall shift below 
T c is naively expected to be A a + A fl ~ 9meV. However, this analysis neglects the multi¬ 
gap nature, which may account for shifts as small as 3meV by allowing scattering into the 
7T band. 

In addition to extracting the energy of the collective mode, we also can estimate the 
electron-phonon coupling \ = iso/is F — l, where u 0 and u F are the bare and renormalized 
Fermi velocity respectively. z/ 0 is estimated from the dispersion at high energy and u F is 
obtained from the dispersion above the kink (as illustrated in Fig. 2a by red dashed and 
blue dotted lines). For cut #2 in Fig. lc, u 0 ~3.77eVA and u F ~1.5 6eVA, which implies 
Ac=1.42. The electron-phonon coupling can be calculated numerically, where it results from 
scattering between the a and 7r-bands, X aw = .23 or from intra-a-band scattering, \ aa = .96, 
for a total of A CT = 1. lf j 19 I 38 I, which is slightly smaller, but within error bars of our results 
given uncertainty of estimating u 0 . 

We now turn to describe the second, low energy excitation present only below the su¬ 
perconducting transition. In Fig. 3 (a-c) we plot the dispersion of the cri-band in close 
proximity to E F below and above T c . At low temperature, well below T c , the shape of the 
dispersion is rather unusual for a superconductor. Instead of back bending of the band as 
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FIG. 2. (Color online) Temperature dependence of the renormalization effect, (a) Band dispersions 
obtained from MDC fits along cut #2 (see Fig. la) measured at different temperatures. Data are 
shifted horizontally for clarity. Band dispersion at 25K (blue solid line) is superimposed with 65K 
data (thick red doted line) for direct comparison of the effect below and above T c . Blue doted and 
red dashed lines illustrate the dispersion used to extract velocity at low and high binding energies 
for estimating coupling constant, (b) Extracted effective real part of self-energy obtained from 
(a). Data are shifted vertically with 25meV interval for each temperature. Red dots mark peak 
position, (c) Temperature dependence of kink energy (red squares) and coupling constant (green 
circles). 
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FIG. 3. (Color online) Low energy band dispersion and kink, (a-c) ARPES intensity measured at 
15K, 29K and 40K along cut #2 (see Fig. la), (d) dispersion at low temperature extracted using 
MDC and EDC fits to data in (a), (e) EDCs measured at 15K, arrow marks location of the dip 
which roughly corresponds to the energy of the Leggett mode, (f) temperature dependence of the 
EDO’s. Clear dip in the spectrum due to interaction with Leggett more is marked by an oval. 


it reaches the energy of the superconducting gap, a sharp streak of intensity is present at 
the gap energy on either side of k^. At 29K (panel b) this feature is also above Ep due 
to increased thermal excitation. All of these features vanish just above T c (panel c), where 
the ordinary conduction band is present. In panel (d) we plot the dispersion of the low 
temperature features extracted using EDC and MDC fits. The sharp peak of intensity is 
almost dispersionless and persists over hk of 0.2 A -1 unlike what is expected for Bogolubov 
quasiparticles 3940 . The EDCs for low temperature data are plotted in panel (e). Here again 
we observe a sharp, dispersion less quasiparticle peak centered at the energy of the SC gap 
and separated from the rest of the spectral weight by a dip at ~10 meV (marked by an 
arrow). The temperature dependence of the EDCs slightly off kp are shown in panel (f). 
A clear dip in the spectrum is observed, which vanishes as the temperature approaches T c . 
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All of these features are due to abrupt changes in the self energy (i. e. onset of resonant 
scattering) and are very characteristic of an interaction between the electrons and a col¬ 
lective excitation. It should be noted that opening of the SC gap can lead to suppression 
of scattering within energy of 3A CT =19.5 meV for intra-band and A a +2A 7T = 9.5 meV for 
inter-band electron channels. We do not observe signatures of suppression of intra-band 
scattering (i. e. line shape features at 19.5 meV) which are deemed to be stronger than 
inter-band ones 41 . Further more, strong suppression of scattering and reduction of E"(c o) 
below certain energy alone cannot produce a spectral dip and non-dispersive peaks. That 
requires a resonant process such as interaction with a collective mode that causes a peak 
in E "(lu). An illustration of this in form of simulations of the spectral function for various 
scenarios of E"(a;) are presented in the Supplementary Material. The lowest energy of an 
optical phonon in MgB 2 is 35 meV and there are no other obvious low energy excitations 
that could couple strongly to the electrons other than the Leggett mode. Whereas Raman 
spectroscopy has found the Leggett mode at 9.4meV^, the similarity of this value to the dip 
energy of 10 meV is coincidental, as the dip occurs at the frequency where electrons above 
the gap may scatter into the collective mode plus quasiparticles above the superconducting 
gap. Therefore we extract the value of the mode energy from ARPES data of 3.5 meV, 
which is the difference between energy location of sharp peak and dip in the spectrum. We 
stress that this is an estimate as the dip location is slightly affected by functional form of 
the self energy. The vertex corrections for Raman spectroscopy 25 30 and ARPES are almost 
certainly different and can potentially explain the the difference between the Leggett mode 
energy measured by the two techniques. Previous electron spectroscopy studies reported 
the value of the Leggett mode energy of 3.9-4.0 meV, consistent with our results. While 
the bare mode was calculated to have a frequency of 5.1-6.2meV 2 ®^, if the experimentally 
measured = 1.5meV is used, the calculated frequencies decrease to 3.8-4.5meV. More 
theoretical efforts will be required to fully understand the origin of the difference. 

We also note that spectral characteristics reported in Fig. 3 (non-dispersive sharp peak 
and presence of a dip) closely resemble features observed at antinode in the cuprates that are 
attributed to strong coupling of electrons to a collective mode. The single order parameter 
of the cuprates does not support a Leggett mode, but these properties are characteristic 
of any collective mode developing only below the superconducing transition. In Fig. 4 we 
show data from optimally doped Bi2212 for comparison. Panels a and b show data at the 
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FIG. 4. (Color online) ARPES data from OP Bi2212 for comparison. Intensity plot at antinode 
in a) normal and b) superconducting state measured using Helium discharge lamp and 21.2 eV 
photons. Dashed line in b) marks location of dip (local minimum of intensity) due to interaction 
with 40 meV collective mode. This feature is absent in normal state (panel a), c) Intensity plot 
and d) EDCs in SC state along cut slightly off the node (marked in inset), where the SC gap 
magnitude is ~8 meV measured using laser source and 6.7 eV photons. The intensity peak above 
E p is due to electrons thermally excited above 2A. 


antinode where electrons couple strongly to a collective mode below T c . In SC state a very 
sharp, weakly dispersing peak emerges at the energy of the SC gap followed by a dip of 
intensity. These two features - lack of substantial dispersion of the sharp quasiparticle peak 
and presence of a dip are key signatures of the interaction with a collective mode. The data 
close to the node (A=8 meV) where the coupling to low energy modes in Bi2212 is very 
weak 9 in SC state display Bogoliubov-like dispersion and sharp, non-dispersive peaks are 
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absent. This is in stark contrast to antinodal direction and present data from MgB 2 . 

Our results demonstrate that in a prototypical, conventional superconductor MgB 2 , the 
coupling of the conduction electrons to phonon mode is not significantly affected by the SC 
transition. This is in contrast to behavior of the collective mode at antinode in cuprates, 
which disappears upon transition to normal state. We also discovered a signature of second 
collective mode in MgB 2 with energy ~ 3.5 meV that exists only below T c . All char¬ 
acteristic of this mode are consistent with Leggett mode which arises due to the relative 
oscillation of the phases of two superconducting condensates present in this material. 
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Supplemental materials for Strong interaction between 
electrons and collective excitations in multiband super¬ 
conductor MgB2 

I. SIMULATION PROCEDURE 

ARPES is an ideal tool with which to examine electron-boson coupling, as both the real, 
E'(cn) and imaginary, E"(A ! ) parts of the self-energy can be straightforwardly extracted from 
momentum distribution curves (MDCs). E ; (u;) is the shift away from the bare dispersion, as 
extrapolated from high energy data, while the scattering rate E"(u;) is proportional to the 
width of the MDC. For interactions with a collective mode at a single frequency, Q, is 

peaked close to Q, while E"(oj) exhibits a step-like increase with the mid point close to Q, 
as quasiparticles can now decay into the mode. The electron-boson coupling strength, A for 
each band depends on the renormalization of the effective mass, A = m*/m — 1. For a more 
complicated phonon spectrum, both parts of the self-energy can simply be calculated from 
the Eliashberg coupling function, a 2 F(u>), which can be calculated numerically, allowing for 
a precise comparision. The temperature dependence of E(u;) should be weak for T <C Q. 
However, as superconductivity gaps out the quasiparticles, the peak/step is expected to 
shift to + A, at least within a simple BCS picture. Superconductivity can similarly af¬ 
fect the energy of the dispersion anomalies, as phonons with energies above (below) the 
gap are expected to harden (soften) when the material is cooled through the transition. 
Strong electron-electron interactions can also strongly renormalize the bandstructure, how¬ 
ever, these effects should be small in MgB 2 . 

Below we illustrate the impact of features in the self energy on ARPES spectra by per¬ 
forming simulations of the ARPES intensity for several scenarios of self energy. We start 
with an assumed form of the imaginary part of the self energy. We then calculate the real 
part of the self energy using a Kramer-Kronig transformation, using a reasonable value of 
the cut-off frequency u> c of leV. The spectral function is then calculated using the following, 
well established formula: 

_ \u k \ 2 E” _ + _ H 2 E” _ 

(a;-S' - y/4 + A 2 ) 2 + E” 2 (w-E+V / 4 + A') 2 + S” 2 
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where k is momentum, u> is binding energy, £' and £” are the real and imaginary parts 
of the self energy respectively, A is the magnitude of the superconducting gap, ek=Vi?k is 
assumed linear band dispersion with Fermi velocity of 4 eVA -1 . Coefficients u/, : and \y. are 
standard BCS coherence factors: 


Kl 2 = 


1 + 






Kl 2 = s i- 


et 


2 v 

The ARPES intensify is then calculated by multiplying the spectral function by the Fermi 
function: J(k, cu) = A(k.,u)f(u). 


II. RESULTS 

A. Normal state 

We begin by plotting in Fig. SI the calculated ARPES intensity and EDCs for the 
normal state assuming the form of self energy based on the data in Fig. 1 of the main text. 
The resulting dispersion is very similar to what is observed in the actual data, with a well 
pronounced kink due the changes in the real part of the self energy. Such changes in the 
real system are usually caused by the onset of scattering. 


B. Superconducting state 

In Fig. S2 we plot the ARPES intensity and dispersion calculated with a 6 meV super¬ 
conducting gap and self energy of a Fermi liquid (£”=a+bn; 2 ). The back bending of the 
band close to the E/ is clearly observed, which is a characteristic signature of Bogolyubov 
quasiparticles. In Fig. S3 we plot the band dispersion and EDCs in the presence of a 6 
meV superconducting gap for a more realistic model of the self energy (similar to Fig. SI 
but with a reduced offset in the imaginary part to better show the features in the proximity 
of E f). The Bogolliubov-like dispersion can be still observed even though the peaks remain 
rather broad close to Ei?. 
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C. Dispersing vs non-dispersing quasiparticle peaks 


The appearance of a sharp quasiparticle peak in the data upon cooling below Tc (Fig. 3 
of the manuscript) signifies a suppression of the scattering rate below a certain energy. We 
attempt to simulate this in Fig. S4, where we include a 6 meV superconducting gap and 
a simple suppression of the scattering rate below 10 meV. This results in a rather abrupt 
sharpening of the quasiparticle peak, but the peaks still display Bogoliubov dispersion, unlike 
what is observed in the experimental data where the sharp peak at the edge of the gap does 
not disperse. This lack of the dispersion requires a ’’resonant” feature in the scattering i. e. 
a peak in the imaginary part, in addition to the suppression of scattering. We illustrate this 
in Fig. S5, which is based on a similar self energy to Fig. S4, but with an additional peak in 
the imaginary part. This peak causes a more rapid change in the real part of the self energy, 
which in turn confines the sharp peak to a small energy range and significantly reduces its 
dispersion. This is of course a grossly exaggerated picture just to demonstrate the effect of 
peaks in the imaginary part of the self energy on the dispersion of the quasiparticles. A 
more realistic model is shown in Fig. S6, where we use the normal state self energy (from 
Fig. SI) that is suppressed below 10 meV and add a small gaussian peak at that energy. 
The resulting APRES intensity and EDO’s very closely resemble the experimental data of 
Fig. 3 in main text. The band dispersion reaches an energy of -10 meV, above which only a 
sharp non-dispersing feature is observed. EDCs show a characteristic sharpening of the line 
shape with the sharp quasiparticle peak sitting at the superconducting gap energy. This 
peak is almost non-dispersing, just like in the data of Fig. 3. The line shape close to kp 
displays a characteristic hump-dip-peak structure, with a dip located close to 10 meV where 
the imaginary part of self energy changes rapidly. 
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FIG. 5. (Color online) Simulated spectral function in the normal state in the presence of a collective 
mode (e. g. phonon) in normal state, a) The imaginary part of the self energy is similar to 
one extracted from the data in Fig. Id of the main text. The real part of the self energy was 
calculated using a Kramers-Kronig transformation, b) The spectral function for the self energy in 
a) is multiplied by the Fermi function, c) same as (b) in close proximity to E f. d) EDCs in the 
proximity of the Fermi wave vector. 


17 















































Energy (meV) 


FIG. 6. (Color online) Simulated spectral function in the superconducting state with a gap of 6 
meV and the self energy of a Fermi Liquid model, a) The imaginary part of the self energy is 
proportional to u 2 with a small offset. The real part of the self energy was calculated using a 
Kramers-Kronig transformation, b) Spectral function for the self energy in a) is multiplied by the 
Fermi function, c) same as (b) but plotted near Ej. d) EDCs in the proximity of the Fermi wave 
vector. Note that the peaks start to develop negative dispersion beyond kp, a key signature of 
Bogolyubov quasiparticle peaks. They lose intensity away from kp due to coherence factors. 
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FIG. 7. (Color online) Simulated spectral function in the superconducting state with a gap of 6 
meV and self energy similar to one in Fig 1, but with a smaller offset, a) The imaginary part 
of the self energy is similar to Fig. la. The real part of the self energy was calculated using a 
Kramers-Kronig transformation, b) Spectral function for the self energy in a) is multiplied by the 
Fermi function, c) same as (b) but plotted near E f. d) EDCs in the proximity of the Fermi wave 
vector. Bogolyubov dispersion is still visible even though the value of ImE close to E^ is large. 
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FIG. 8. (Color online) Simulated spectral function in the superconducting state with a gap of 
6 meV and self energy that is strongly suppressed below 10 meV. a) Imaginary part of the self 
energy. The real part of the self energy was calculated using a Kramers-Kronig transformation, b) 
Spectral function for the self energy in a) is multiplied by the Fermi function, c) same as (b) but 
plotted near E f. d) EDCs in the proximity of the Fermi wave vector. Bogolyubov dispersion is 
clearly observed above 6meV. 
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FIG. 9. (Color online) Simulated spectral function in the superconducting state with a gap of 6 
meV and exaggerated self energy that is strongly suppressed below 10 meV with an additional large 
gaussian peak at 10 meV. a) Imaginary part of the self energy. The real part of the self energy was 
calculated using a Kramers-Kronig transformation, b) Spectral function for the self energy in a) is 
multiplied by the Fermi function, c) same as (b) but plotted near E j. d) EDCs in the proximity 
of the Fermi wave vector. Bogoliubov dispersion is not observed, as the rapidly changing real part 
of the self energy confines the sharp quasiparticle peak. 
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FIG. 10. (Color online) Simulated spectral function in the superconducting state with a gap of 6 
meV and a realistic self energy based on data from Fig. Id of the main text, but suppressed below 
10 meV and with a small gaussian peak at that energy, a) Imaginary part of the self energy. The 
real part of the self energy was calculated using a Kramers-Kronig transformation, b) Spectral 
function for the self energy in a) multiplied by the Fermi function, c) same as (b) but plotted near 
E f. d) EDCs in the proximity of the Fermi wave vector. This simulation looks very similar to data 
in Fig. 3. The sharp peak present at the superconducting gap value is almost non-dispersing. The 
lineshape has a characteristic hump-peak-dip structure. 
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